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Abstract 

Knowledge of both the mass and radius of an exoplanet allows us to estimate its 
mean density, and therefore its composition. Exoplanets seem to fill a very large pa¬ 
rameter space in terms of mass and composition, and unlike the solar-system’s planets, 
exoplanets also have intermediate masses (~5 - 50 M®) with various densities. In 
this letter, we investigate the behavior of the Mass-Radius relation for methane (CH 4 ) 
planets and show that when methane planets are massive enough (M p > 15 M®), 
the methane can dissociate and lead to a differentiated planet with a carbon core, a 
methane envelope, and a hydrogen atmosphere. The contribution of a rocky core to the 
behavior of CH 4 planet is considered as well. We also develop interior models for several 
detected intermediate-mass planets that could, in principle, be methane/methane-rich 
planets. The example of methane planets emphasizes the complexity of the Mass- 
Radius relation and the challenge involved in uniquely inferring the planetary compo¬ 
sition. 


1 Introduction 

Exoplanet studies are now reaching the level at which planet characterization is possible. 
There are hundreds of planets with measured mass and radius, not including the thousands 
of candidates from the Kepler mission. Characterization of exoplanets is often performed via 
the Mass-Radius relationship (hereafter, M-R relation) which provides information on the 
planet’s mean density, and therefore its composition. The curves of planetary radius versus 
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mass for various compositions have been computed by several groups (e.g., Valencia et al., 
2006; Sotin et ah, 2007; Seager et al. 2007; Fortney et al., 2007; Rogers et al., 2011; Lopez & 
Fortney 2014; Lissauer et al., 2011; Fressin et al., 2012; Wagner et al., 2011). The increasing 
number of observed exoplanets with masses of 5 - 20 M e and the large variation in their 
mean densities encourages the investigation of compositions that include volatile materials 
and various mixtures. 

Typically, the M-R curves are smooth and monotonic. Discontinuities in the M-R curve 
could, however, appear when phase changes are considered. An example of this situation can 
occur in hydrogen-dominated planets, where at pressure of ~100 GPa molecular hydrogen 
transitions to metallic hydrogen. This transition leads to a kink in the M-R curve for planets 
with masses of 100 M e (e.g., Swift et ah, 2012). 

Discontinuities in the M-R relation can also occur for other materials. One example is 
methane (CH4). Here, we show that CH4 planets can dissociate and lead to a formation of a 
carbon (diamond) core, a CH4 envelope, and an H2 atmosphere. Under such a configuration, 
the planetary radius can be significantly larger than that expected without dissociation, 
which can affect the composition inferred from the M-R relation. CH 4 dissociation leading 
to a formation of a carbon core has already been presented in the context of the solar-system 
giant planets (Ross, 1981; Benedetti et ah, 1999). Similar behavior can occur for other 
molecules, but in this paper we investigate only CH 4 . Methane is a good example because 
this molecule is expected to be abundant in planetary systems. In the solar-system, for 
example, all of the outer planets are found to be enriched with carbon (by measuring the 
carrier CH 4 ) with Jupiter and Saturn being enriched with carbon by factors of about 4.3 
and 9.7 compared to the proto-solar value, respectively, and Uranus and Neptune by factors 
of about 86 and 68, respectively (e.g., Guillot & Gautier, 2014), although these numbers 
correspond only to the atmospheric composition and it is unclear whether they represent 
the bulk compositions of the planets. Clearly, other molecules such as H 2 0, NH 3 , and PH 3 
could exist as well, and in fact, a planet that is enriched with carbon is likely to have much 
oxygen as well. However, modeling CH4 planets is already non-trivial due to the possibility 
of dissociation and differentiation. 
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Finally, while the behavior of refractory materials such as silicon and iron is relatively 
insensitive to temperature, this is not the case for volatile materials for which could increase 
the planetary radius considerably with increasing temperature. Below, we explore how the 
choice of the dissociation pressure and internal temperature affect the M-R relation for 
methane planets, and investigate how the inclusion of a silicate (or iron) core changes the 
results. Finally, we identify several detected exoplanets that could be methane-rich planets. 

2 Modeling Methane Planets 

2.1 Computing the internal structure 

The internal structure is given by the basic equations of planetary structure: the equation 
of hydrostatic equilibrium and the equation of mass conservation. For simplicity, an energy 
transport equation is not included and the planets are assumed to be isothermal. 

The structure equations are integrated using a fourth-order Runge-Kutta procedure. We 
integrate from a central pressure P c to a low surface pressure P s (typically set to 1 bar) and 
compute the planetary mass and radius. For pure CH 4 planets, if P c is greater than the 
methane dissociation pressure, hereafter, Pdiss, we use the equation of state (EOS) for pure 
carbon. The integration is continued until the pressure falls below Pdiss- The mass of the 
body at this point is the mass of the carbon core, Me- The dissociation of CH 4 produces 
an additional 1/3 Me of hydrogen which we assume has diffused to the outer envelope of the 
planet to form a hydrogen atmosphere. Any remaining mass is in the form of an intermediate 
shell of CH 4 . 

In principle, each planetary mass should correspond to a unique P c value, and we can 
find a single-valued function M (P c ) for the planetary mass versus central pressure. However 
this is not always the case. In the immediate vicinity of P c = Pdiss, central pressures just 
below and just above Pdi SS give the same mass (e.g., Ramsey, 1948). In this case, the planet 
with a carbon core is the stable configuration. 

For planets with a silicate/iron core (see section 2.4) we follow the same procedure, except 
that the innermost region is assumed to be SiC^/Fe. We integrate out from the center until 
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the mass of the body equals the mass of the prescribed silicate/iron core. After that, if the 
pressure is greater than Pdiss we use the equation of state for carbon, otherwise we use the 
EOS for CH4. Also here we guarantee that the mass of the hydrogen atmosphere is one-third 
the mass of the carbon shell. 

2.2 The equation of state 

We use the SESAME EOS for carbon, methane, and hydrogen (e.g., Kerley, 1972). When 
cores are included, we use the SESAME EOS for Si 02 and iron (Fe). In order to account 
for the effect of methane dissociation, we must assume the pressure at which dissociation 
occurs. Several studies suggest that methane dissociates at a pressure of ~300 GPa (Chau 
et al., 2011 ; Gao et ah, 2010 ; Ancilotto et ah, 1997), thus, the exact dissociation pressure 
of CH 4 is still unknown and needs to be confirmed through experiments and theoretical 
calculations. In order to investigate the sensitivity of the derived internal structure (and 
M-R relation) to the dissociation pressure, we also consider a lower dissociation pressure of 
170 GPa. The lower the dissociation pressure is, the lower the planetary mass for which a 
jump occurs in the M-R relation. In fact, before CH 4 breaks down to carbon and hydrogen 
some intermediate products of hydrocarbons can be formed: methane will Erst form ethane 
(C 2 H 6 ) at 95 GPa, butane (C 4 H 10 ) at 158 GPa, and finally carbon (diamond) and hydrogen 
(Gao et ah, 2010). 

If hydrogen diffuses through the shell of undissociated CH 4 to the planetary atmosphere, 
then pure CH4 planets above a certain mass would have a pure carbon core, a methane inner 
shell, and a (molecular) hydrogen atmosphere. In this work, we assume that once dissociation 
occurs differentiation takes place. Modeling the actual of differentiation is beyond the scope 
of this paper. Some pioneering work in that direction was recently presented by Levi et 
al. (2014). 

2.3 Pure methane planets 

For a pure methane planet, there is a sharp discontinuity in the M-R relation above the 
mass where dissociation occurs. This is presented in the two top panels of Figure 1. Since 
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unlike for refractory materials, the temperature is not negligible, we consider four different 
temperatures: 50 K, 500 K, 1,000 K, and 1,500 K. The isotherm of 50K should be taken 
as extreme-low-temperature while other isotherms are comparable to the atmospheric tem¬ 
peratures of close-in exoplanets. Of course, planets are not isothermal but are expected to 
have a temperature profile, with the temperature increasing toward the center. Thus, since 
the temperature mostly affects the planetary radius when dissociation occurs due to the 
formation of an hydrogen envelope, assuming isothermal structure is reasonable. The effect 
of temperatures is much less profound in the deep interior that is composed of refractory 
materials. In addition, the atmospheres of close-in planets are expected to be nearly isother¬ 
mal down to a depth of ~ 100 bar (e.g., Madhusudhan & Seager, 2010), the bottom of the 
hydrogen envelope is convective, in this case an adiabat is more appropriate but since the 
region between 100 bar and the hydrogen-methane boundary is relatively small an adiabat 
profile for the temperature only increases the temperature by a factor of ~ three. In fact, 
since in our case the atmosphere is composed of pure hydrogen the opacity is expected to be 
lower than the typically modeled atmospheres since they include H 2 0, CO, CH 4 and C0 2 
that increase the opacity, and therefore the radiative-convective boundary of pure-hydrogen 
atmospheres, as in our case, will be at even higher pressures. We can then conclude that an 
isothermal structure is a reasonable assumption but clearly the sensitivity of the assumed 
temperature profile on the structure should be investigated in the future. The effect of tem¬ 
peratures on the M-R relation is strongest in the outer region (hydrogen) and is much less 
profound in the deep interior that is composed of refractory materials. 

The derived M-R relations are shown for P diss of 170 GPa (left) and 300 GPa (right). The 
solid, dashed, dotted, and dashed-dotted curves correspond to T = 50 K, 500 K, 1,000 K 
and 1,500 K, respectively. As long as the planet is undifferentiated, the radius is insensitive 
to temperature. When the planetary mass is high enough to produce a hydrogen (H 2 ) 
atmosphere, the radius increases abruptly and becomes much more sensitive to temperature. 
The amount of methane in the planet varies with its mass. As long as the planetary mass 
is low enough so that dissociation does not occur, the mass fraction of methane is unity 
and the radius is insensitive to the temperature, thus, once the critical mass (pressure) for 
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dissociation is reached, the mass fraction of methane decreases rapidly. 

Hot methane planets that dissociate and develop a molecular hydrogen envelope, could, 
in principle, lose some of their atmospheres. For a sun-like star, the four temperatures of 
50 K, 500 K, 1,000 K and 1,500 K correspond to radial distances of about 30 AU, 0.33AU, 
0.083 AU and 0.037 AU, respectively. The most vulnerable planets are those with the high 
temperatures (i.e., T = 1,000 K and T = 1,500 K). Planets with small masses have P c < Pdiss 
and therefore are not expected to have hydrogen atmospheres. For the cases we investigate, 
this means masses less than ~ 8 M e (see Fig. 1). More massive planets can dissociate CH 4 , 
but for these masses hydrogen escape is less efficient due to the relatively large gravity. We 
can estimate the mass-loss rate for the different cases by using the energy-limited mass loss 
estimates (e.g., LeCavelier des Etangs, 2007; Lopez et al., 2012, Owen & Jackson 2012; 
Murry-Clay et al., 2009; Mordasini et ah, 2015; Kurosaki et ah, 2014). For planets with 
temperatures of 50 and 500 K evaporation is found to be negligible. For temperatures of 
1,000 K and 1,500 K, planets with masses larger than ~ 7 M® the mass that is lost after 
10 Gyr is about 2 % and 10 % of the total planetary mass, respectivel}Q Smaller planets, 
however, can lose a significant amount of mass, especially if they are close to their host star. 
Therefore, low-mass planets with hydrogen atmospheres are unlikely to exist in old systems. 
Thus, the mass loss estimate assumes a pure-hydrogen atmosphere, and since small planets 
cannot dissociate (due to low internal pressures), they have methane atmospheres, for which 
the mass loss rate can be much lower. Further investigations of the mass loss rates of various 
atmospheric compositions is required to understand the survivability of such planets (e.g., 
Lammer et al., 2012). 

2.4 The effect of a silicate core 

Since planets are not expected to consists of pure methane, we next consider the existence of 
a silicate (SiCU) core and investigate how it affects the derived M-R relation. The presence 

lr The estimates were performed assuming energy-limited mass loss around a G star, and accounting for 
the change in the X-ray and UV radiation flux at stellar age of 0.1 Gyr (Ribas et al., 2005) integrated up 
to an age of 10 Gyr. More information can be found in Kurosaki et al. (2014) and references therein. These 
estimates however are very crude since the change in internal structure due to evaporation is not included. 
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of a silicate core reduces the importance of this effect. As the silicate mass fraction of the 
planet increases, the magnitude of the radius change decreases. This is because the silicate 
core occupies part of the high pressure region where the methane dissociation should take 
place. The larger the silicate core, the smaller the volume of methane that participates in 
the dissociation, and the smaller the difference in radius at the transition point. The effect 
of the silicate core is demonstrated in the bottom panels of Figure 1 where we compare the 
M-R relation for methane planets with no core (solid-black) and with silicate cores that are 
equal to 30% (dotted-gray), 50% (dashed-blue) and 80% (dashed dotted-red) of the total 
planetary mass for T= 500 K for Pdiss of 170 GPa (left) and 300 GPa (right). While we 
use SiC >2 to represent the core material, the core could be composed of other elements such 
as Fe or MgO. The results are not expected to change much qualitatively since the effect is 
similar as long as the core material is significantly denser than methane. A comparison of 
the inferred internal structure for an SiC >2 core vs. an Fe core is discussed in the following 
section. As can be seen from the bottom panels of Figure 1, the discontinuity due to CH 4 
dissociation is much less prominent when the planet consists of a massive core. When the 
silicate core consists of 30% and 50% of the total planetary mass, dissociation occurs at 
smaller masses. This is because the core results in a higher internal pressure which allows 
dissociation at lower masses than in the no-core case. Thus, at some point, when the core 
mass is significantly larger than the methane mass, methane will exist only at the outer 
region (envelope) where there pressure is low. In that case, dissociation occurs at higher 
masses compared to the no-core case. Indeed, for the case of 80% core dissociation is found 
to occur at higher masses. Here too the kink in the M-R relation occurs for larger masses 
for the higher Pdiss- 

3 Fits to observed exoplanets 

Several detected exoplanets have relatively low masses but radii that are larger than expected 
from pure H 2 0, and therefore must consist of an atmosphere. Here, we show that some of 
these planets are good candidates for methane-dominated planets. In Figure 2 we show the 
M-R relation for planets with masses up to 20 M e for various materials including several 
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detected exoplanets. The curves for pure iron (brown), 100% MgSi0 3 (black), H 2 0-MgSi0 3 
(green), and 100% water (cyan) are taken from Zeng & Sasselov (2013). Also shown are the 
M-R curves for methane/methane-rich planets. Several of the detected exoplanets have, in 
principle, masses and radii that are consistent with those predicted for methane/methane- 
rich planets. Clearly, there are other possible solutions and as argued in this Letter, there 
is no unique solution even for methane planets. Nevertheless, the possibility of methane- 
dominated planets is intriguing. 

Examples for potential-methane planets are GJ3470b and Kepler 89c[^] The inferred 
internal structure for these planets are shown in Figure 3 (left panel). The inferred interior 
structure is shown for four different temperatures: 50 K (blue), 500 K (gray) 1,000 K (black), 
and 1,500 K (red). For colder planets, the mass fraction of the silicate core is larger, but in all 
the cases the planet has a differentiated structure. For Kepler-89c, the region of the methane 
shell is found to be significantly larger for Pdiss = 300 GPa, and the hydrogen atmosphere 
is more compact. In order to demonstrate the sensitivity of the internal inferred structure 
to the core material we compare the solutions for T = 1,000 K for Si0 2 and iron (Fe) cores 
(right panels). 

The derived values for the mass fractions and locations for the different elements are 
summarized in Table 1, also listed are other exoplanets for which there is a solution consistent 
with differentiated methane planets. Examples for exoplanets that are consistent with an 
internal structure which consists of Si0 2 and CH 4 but with undissociated methane are listed 
in Table 2. The uncertainty in the measured radii and masses can be very large. Here, 
we set the mass and radius of a given exoplanet to be their 0-sigma values. Clearly, the 

2 The estimate for the mass and radius used in this paper for Kepler-89c were based on Weiss et al. (20f 3) 
with M = 15.57i5 5 7 f M® and R = 4.33 ± 0.09 R®. New estimates for the mass and radius of the planet 
were derived by Masuda et al. (2013), who suggest a planetary mass of 9.lb 2 '] M® and radius of 3.8b2 6 29 R®- 
If Kepler-89c is indeed less massive and smaller as suggested by Masuda et al. (2013), Kepler-89c will not be 
able to differentiate for a dissociation pressure of 300 GPa, and in addition, for the cases were differentiation 
does occur (170 GPa) the internal structure of the planet will change as well. For the numbers suggested 
by Masuda et al. (2013) a smaller fraction of methane can dissociate, and therefore, Kepler-89c will have 
smaller carbon and molecular hydrogen shells than presented here. Thus, Kepler-89c still remains a potential 
differentiated-methane-rich planet for Rjiss = 170 GPa. 



uncertainty in mass and radius introduce a large uncertainty on the actual internal structure 
of the planet. This emphasizes the need for accurate determinations of the masses and radii 
of exoplanets. 

4 Discussion 

Observations of exoplanets have taught us that planets come in a variety of masses and com¬ 
positions. While the M-R relation is used for classifying exoplanets, the relation between 
mass and radius is not always straightforward. Here, we show that even a simple compo¬ 
sition consisting of one species, such as methane, can lead to complications due to phase 
transition and differentiation. In addition, the temperature must be accounted for when the 
exoplanets consist of volatile materials. We show that some of the observed masses and radii 
of intermediate-mass exoplanets are, in principle, consistent with methane-rich interiors. 

Our model is rather simple. A necessary next step would be to derive M-R relations 
for planets with more realistic temperature profiles (e.g., an adiabatic temperature profile). 
In addition, future models should account for other possible mixtures. In particular, it is 
interesting to investigate how the presence of oxygen affects the derived internal profile, 
and to investigate whether or not the existence of other elements prevents the formation of 
a hydrogen atmosphere. The presented research can be continued by accounting for more 
materials, but for that case chemical interactions should be modeled as well. Finally, a 
better theoretical understanding of the differentiation of methane planets and knowledge 
of the behavior of methane at high pressure and temperatures are required. Upcoming 
space missions such as TESS, CHEOPS and PLATO 2.0 and ground-based observations will 
provide new and exciting data, and it is now an ideal time to establish a deeper understanding 
of the M-R relation and its potential for characterizing exoplanets. 
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Figure 1: Top: M-R relation for pure CH 4 planets with isotherms of 50 K (solid), 500 K (dotted), and 
1,000 K (dashed), and 1,500 K (dashed-dotted) assuming a dissociation pressure of 170 GPa (left panel) and 
300 GPa (right panel). Bottom: M-R relation for methane planets with different mass fractions of a silicate 
(Si 02 ) core for a temperature of 500 K. Also here the left and right panels correspond to Pdiss of 170 GPa 
and 300 GPa, respectively. See text for details. 


12 

















Mass (M m ) 

Figure 2: M-R relation for various compositions including methane, and several exoplanets with detected 
masses and radii. The curves for H 2 0, Fe-MgSiC> 3 , MgSiC >3 are taken from Zeng & Sasselov (2013). The 
black curves correspond to pure CH 4 using Pdi ss = 170 GPa for the four temperatures we consider (similar 
to Fig. 1, top-left). The gray curves show the M-R curves for methane planets with silicate (Si0 2 ) cores 
consisting of 30% (dotted) and 50% (dashed) of the total planetary mass. Several planets are found to have 
masses and radii that are consistent with pure methane planets, such as Kepler 89c, GJ3470b, and three of 
the Kepler 11 planets: Kepler 11c, Kepler lid, and Kepler Ilf. The masses and radii of Kepler lib and GJ 
1214b fit a mixture of methane and rocks. 
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Figure 3: Left: Density vs. normalized radius for GJ3470 b with P d i ss =170 GPa (top), Kepler 89c with 
Pdiss=170 GPa (middle), and Kepler 89c with P diss =300 GPa (bottom). The solid-black, dotted-gray, dashed- 
blue, and dashed-dotted-red curves correspond to temperatures of 50 K, 500 K, 1,000 K, and 1,500 K, 
respectively. Right: Density vs. normalized radius for GJ3470 b with P diss =170 GPa (top), Kepler 89c with 
Pdiss=170 GPa (middle) and Kepler 89c with P diss =300 GPa (bottom) for T = 1,000 K for two different 
core materials: Si02 (solid-brown) and Fe (dotted-black). As can be seen from the figure, an iron core leads 
to much higher internal densities. 
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Table 1: Solutions for a few detected exoplanets that are consistent with differentiated methane planets. 
Msio 2 /M p , Mc/Mp, Mch 4 /M p , Mh 2 /M p correspond to the mass fractions of Si0 2 , C, CH4 and H 2 , respec¬ 
tively. Rsi 0 2 /Rp, Rc/R-p, RcH 4 /Rp, and Rh 2 /R p correspond to the radii (normalized to the planet’s radius) 
of the shells for the compositions of Si0 2 , C, CH 4 and H 2 , respectively. Given are the solutions for the 
calculated interior structure assuming dissociation pressures of 170 and 300 GPa, and four different internal 
temperatures: 50 K, 500 K, 1,000 K, and 1,500 K. The masses and radii of the planets are taken from the 
exoplanet database at http://exoplanets.org and are taken to be their Ocr values. 
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Undifferentiated 

Kepler lid 

M=6.1 M e 

R=3.117 ±°;“ R e 




Msio 2 /Mp 

RsiC^ /-^p 

Mch 4 /M p 

Rch 4 /R p 

50 K 

0.18 

0.354 

0.82 

1.0 

500 K 

0.22 

0.376 

0.78 

1.0 

1,000 K 

0.28 

0.403 

0.72 

1.0 

1,500 K 

0.33 

0.423 

0.67 

1.0 

GJ 1214b 

M=6.47 ± 1.0 M e 

R=2.68 ± 0.13 R 0 




Msio 2 /Mp 

R-Si02 

Mch 4 /M p 

Rch 4 /R p 

50 K 

0.50 

0.544 

0.5 

1.0 

500 K 

0.53 

0.554 

0.47 

1.0 

1,000 K 

0.58 

0.570 

0.42 

1.0 

1,500 K 

0.62 

0.583 

0.38 

1.0 

HD 97658b 

M=7.87 ± 0.73 M e 

R=2.34 t° 0 il R® 




Msio 2 /Mp 

-f^Si02 /-^p 

McH 4 /Mp 

Rch 4 /Rp 

50 K 

0.88 

0.811 

0.12 

1.0 

500 K 

0.895 

0.818 

0.105 

1.0 

1,000 K 

0.91 

0.824 

0.09 

1.0 

1,500 K 

0.92 

0.83 

0.08 

1.0 

Kepler 68b 

M=6.0 ± 1.71 M e 

R=2.308 t°o°ot R® 




Msio 2 /Mp 

Rsi02 /f^p 

Mch 4 /M p 

Rch 4 /Rp 

50 K 

0.85 

0.782 

0.15 

1.0 

500 K 

0.865 

0.788 

0.135 

1.0 

1,000 K 

0.88 

0.795 

0.12 

1.0 

1,500 K 

0.89 

0.799 

0.11 

1.0 


Table 2: Solutions for exoplanets that potentially have methane-rich interiors but for which dissociation of 
methane does not occur. Shown are the normalized masses and radii of the silicate (SiOo) core and methane 
(CH 4 ) envelop for the four different internal temperatures used in this work. Msio 2 /Mp and Mch 4 /M p 
correspond to the mass fractions of SiC >2 and CH 4 , respectively, while Rsi 0 2 /Rp andRcH 4 /Rp are the ra dh 
(normalized to the planet’s radius) of the shells for the compositions of SiC> 2 , and CH 4 , respectively. Also 
here, the masses and radii of the planets are taken from the exoplanet database at http://exoplanets.org 
and are taken to be their Ocr values. 
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